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1.  IntroductKMi 


This  is  the  fined  report  on  the  URI  project  entitled  “Computational  £uid 
Analytical  Methods  in  Nonlinear  Fluid  Dynamics”  started  in  1989  at  Lehigh  University 
in  the  Department  of  Mechanical  Engineering  and  Mechanics.  The  present  report 
summarizes  work  completed,  as  well  as  still  in  progress,  at  the  end  of  the  funding.  The 
central  focus  of  the  program  was  on  the  application  and  development  of  modern 
analytical  and  computational  methods  to  the  solution  of  nonlinear  problems  in  fluid 
dyn2unics  and  reactive  gas  dynamics.  The  research  was  carried  out  within  the  Division 
of  Engineering  Mathematics  in  the  Department  of  Mechanical  Engineering  and 
Mechanics  and  principally  involved  Professors  P.  A.  Blythe,  E.  Varley  and  J.  D.  A. 
Walker.  In  addition,  the  program  involved  various  international  collaborations. 
Professor  Blythe  completed  work  on  reactive  gas  dynamics  with  Professor  D.  Crighton 
F.R.S.  of  Cambridge  University  in  the  United  Kingdom.  Professor  Walker  and  his 
students  carried  out  joint  work  with  Professor  F.  T.  Smith,  F.R.S.  of  University  College 
London  on  vcirious  problems  in  unsteady  flow  and  turbulent  boundary  layers.  This 
research  will  be  briefly  summarized  subsequently. 

A  major  portion  of  the  work  carried  out  under  the  contract  was  computational 
and  concerned  the  incorporation  of  analytical  local  solutions  into  computational 
algorithms.  In  a<.cordance  with  the  theme  of  the  URI  program  to  enhance  University 
infrastructure,  two  Silicon  Graphics  4D25  workstations  were  purchased  with  contract 
funds.  These  computers  were  connected  to  the  University  network  in  a  computational 
lab  principally  assigned  to  fluid  mechanics  research.  They  received  extensive  usage  over 
the  course  of  the  contract  and  are  still  regularly  used  on  a  daily  basis  by  graduate 
students  and  faculty. 

During  the  contract  period,  four  graduate  students  were  funded  (in  full  or  in 
part)  for  studies  toward  their  Ph.D.  These  students  were  supervised  by  Professor 
J.  D.  A.  Walker  and  at  present  one  Ph.D.  has  been  awarded;  the  three  other  students 
will  receive  their  Ph.D’s  in  December,  1993.  The  students  and  their  thesis  topics  are  as 
follows: 

(1)  Ammar  T.  Degani,  Ph.D.  1991,  “The  Three-Dimensional  Turbulent  Boundary 
Layer  -  Theory  and  Application”. 

(2)  Kevin  Cassel,  Ph.D.  (expected)  1993,  “Interactive  Boundary  Layers”. 
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(3)  Jun  He,  Ph.D.  (expected)  1993,  “High-Speed  Compressible  Turbulent  Boundary 
Layers” . 

(4)  R.  I.  Puhak,  Ph.D.  (expected)  1993,  “Flow  Structure  Upstream  of  Three- 
Dimensional  Obstacles”. 


A 


Mr.  Cassel  and  Mr.  Puhak  are  American  citizens  while  Dr.  Degani  is  a 
residential  alien.  Dr.  Degani  recently  won  a  Nason  Postdoctoral  Fellowship  in  Parallel 
Processing  at  Syracuse  University,  while  Mr.  Cassel  was  successful  in  a  competition  for 
a  Postdoctoral  Fellowship  at  NIST  in  Combustion  Studies. 

In  the  following  sections,  research  carried  out  under  the  gr<mt  is  briefly 
summarized;  more  detail  is  available  in  the  publications  summarized  at  the  end  of  this 
report. 

2.  Description  of  Research 

A.  Professor  Di  Ai  Walker 

A.l  Analytical  and  Computational  Studies  of  Unsteady  Separation 
Associated  Graduate  Student:  Kevin  Cassel 

Status:  Ph.D.,  draft  of  thesis  complete,  expected  graduation  December  1993. 

Unsteady  separation  of  boundciry  layers  occurs  in  a  variety  of  physical  situations  • 

such  2is  dynamic  stall  in  rapidly  pitching  airfoils,  vortex  shedding  in  turbines  and 
eruption  of  the  wall  layer  in  transitioning  and  turbulent  boundary  layers.  Until 
relatively  recently,  the  nature  and  causes  of  this  phenomena  at  high  Reynolds  numbers 
have  not  been  well  understood.  The  Ph.D.  thesis  of  Van  Dommelen  (1981)  and  • 

subsequent  theoretical  work  by  Elliott  et  al.  (1983)  revealed  that  in  two-dimensional 
flows,  a  generic  response  develops  within  a  thin  boundary- layer  flow  that  is  exposed  to  a 
persistent  adverse  pressure  gradient.  The  influence  of  viscosity  first  leads  to  the 
evolution  of  recirculating  flow  within  the  boundary  layer  and  if  the  external  julverse  • 

pressure  gradient  is  maintained,  a  sharply  focussed  eruption  of  the  surftice  layer 
develops  abruptly  at  a  subsequent  time.  The  calculation  of  such  flows,  using 
conventional  numerical  methods  based  on  the  Eulerian  description  of  the  flow  field,  is 
not  possible  owing  to  the  rapid  development  of  a  “spike-like”  eruption  of  the  surface  • 

flow.  Van  Dommelen  (1981)  was  able  to  calculate  the  evolution  of  the  boundary  layer 
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on  an  impulsively-moved  circular  cylinder  through  the  use  of  Lagrangian  methods, 
wherein  the  trajectories  of  a  large  number  of  individual  fluid  particles  were  calculated. 
In  this  manner,  it  proved  possible  to  track  the  phenomenon  up  to  the  formation  of  a 
sharp  spike  in  displacement  thickness.  The  theoretical  account  of  Elliott  et  al.  (1983) 
suggests  that  the  formation  of  the  “spikes”  is  a  generic  state  passed  through  by  most 
erupting  boundary  layers  at  high  Reynolds  numbers.  Recently,  this  suggestion  has  been 
supported  by  numerical  solutions  of  the  boundary  layer  induced  by  a  vortex  above  a 
plane  wall,  carried  out  by  Peridier  et  al.  (1991a,  1991b).  In  the  limit  of  large  Reynolds 
numbers,  the  vortex-induced  separating  boundary  layer  evolves  toward  a  singularity 
structure  which  is  essentially  similar  to  that  described  by  Elliott  et  al.  (1983).  As  the 
singularity  structure  develops,  the  external  flow  field  must  be  allowed  to  change  and 
reflect  the  fact  that  boundary- layer  vorticity  is  now  shooting  rapidly  toward  the  outer 
region. 

Because  unsteady  separation  at  high  Reynolds  numbers  is  pervasive  in 
aeronautical  applications,  it  is  important  to  develop  an  understanding  of  the  steps  that 
occur  in  the  process  as  a  separating  boundary  layer  leaves  the  surface.  Until  recently,  it 
was  believed  by  many  authors  that  the  problem  could  be  haindled  using  conventional 
interacting  bound^y  layer  concepts.  The  present  research  has  definitively  shown  that 
this  is  not  the  case  and  that  more  sophisticated  algorithms  will  be  required  for  this 
importcint  problem.  In  1983,  Elliott  et  al.  formulated  a  boundary  value  problem 
describing  what  we  refer  to  as  the  first  interactive  stage:  this  occurs  just  prior  to  the 
evolution  of  a  singularity  in  the  non-interactive  problem  and  is  the  first  stage  wherein 
the  external  flow  is  influenced  by  the  rapidly  thickening  viscous  region  near  the  surface. 
The  unsteady  boundary- value  associated  with  this  stage  is  extremely  complex  from  a 
numerical  standpoint  owing  to  difficulties  associated  with  the  boundary  conditions. 

The  first  numerical  solutions  for  the  first  interactive  stage  were  obtained  in  the 
present  program.  By  rec^lsting  the  problem  in  a  Lagrangian  formulation,  it  proved 
jx)ssible  to  formulate  a  numerical  cilgorithm  that  would  accurately  calculate  the  flow 
development  in  the  first  interactive  stage.  Initially,  it  was  expected  that  the  solution 
would  evolve  toward  another  singularity  but  repeated  calculations  revealed  the 
evolution  of  a  high  frequency  instability.  As  the  meshes  were  refined,  this  instability 
was  found  to  occur  at  progressively  earlier  stages.  Repeated  checks  on  the  numerical 
work  and  analysis  of  the  instability  were  carried  out.  This  work  has  shown  that  the 
instability  is  not  numerical  but  appears  to  be  a  new  type  of  interactive  instability.  It  is 
primarily  inviscid  in  nature  but  does  not  involve  an  inflectional  velocity  profile.  A 
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paper  on  this  discovery  has  been  submitted  to  the  Journal  of  Fluid  Mechanics  (Cassel, 
Smith  &  Walker,  1993). 

A.2  Embedded  Function  Methods  for  Three-Dimensional  Turbulent  Boundary  Lavers 
Associated  Grciduate  Student:  A.  T.  Degani 
Status:  Ph.D.  awarded  in  October,  1991 

In  most  conventional  methods  for  the  calculation  of  three-dimensional  turbulent 
boundary  layers,  the  majority  of  computational  time  and  resources  is  invested  in 
resolving  the  severe  gradients  in  the  velocity  and  temperature  gradients  that  occur  in 
the  near-wall  region.  In  the  present  research  asymptotic  analysis  w<is  used  to  develop 
local  ein£dytical  solutions  that  could  be  embedded  in  a  general  computational  scheme  for 
the  turbulent  boundary  layer.  In  the  papers  by  Degani  et  al.  (1992,  1993),  the  method 
of  macched  expansions  was  used  to  determine  the  general  asymptotic  structure  of  an 
attached  three-dimensional  turbulent  boundary  layer.  This  work  established  for  the 
first  time  that  the  inner  wall  layer  exhibits  an  essentially  universcd  self-similar  behavior 
and  the  matching  conditions  for  the  velocity  profiles  near  the  edge  of  the  wall  layer 
were  determined.  A  set  of  analytical  profiles  that  can  be  used  in  the  turbulent  wall 
layer  in  a  computational  algorithm  were  developed.  The  concept  was  demonstrated  in  a 
three-dimensioned  boundary- layer  prediction  code  (Degani  &  Walker,  1993)  where  the 
analytic  embedded  functions  were  used  to  “step  over”  the  turbulent  wall  layer  in  a 
calculation.  In  this  manner,  it  is  only  necessary  to  carry  out  a  computation  in  the  outer 
part  of  the  boundary  layer  and,  skin  friction  and  heat  transfer  coefficients  are 
determined  through  matching  to  the  analytic  embedded  functions.  Through  use  of 
asymptotic  theory  and  the  embedded  function  methodology,  it  h^ls  proved  possible  to 
construct  a  much  more  efficient  algorithm  for  turbulence  prediction  near  w6Jls  in  which 
the  total  number  of  mesh  points  used  may  be  reduced  by  about  50%  (as  compared  to  a 
conventional  method)  with  no  degradation  in  accuracy. 

A. 3  Modeling  of  Supersonic  and  Hypersonic  Turbulent  Boundary  Lavers 
Associated  Graduate  Student:  Jun  He 

Status:  Thesis  draft  complete,  Ph.D.  expected  December,  1993. 

Supersonic  and  hypersonic  turbulent  boundary  layers  occur  in  a  variety  of 
import2mt  aerodynamic  applications.  Unfortunately,  most  turbulence  models  for  high¬ 
speed  compressible  flow  are  based  upon  ad  hoc  extensions  of  turbulence  models  that 
have  been  correlated  to  represent  subsonic  data.  As  the  Mach  number  of  the 
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mainstream  flow  increases,  the  structure  of  the  turbulent  boundary  changes  and  a  new 
physical  effect  becomes  important,  namely  dissipation  of  kinetic  energy  near  the  surface 
into  heat  energy.  Unfortunately,  there  are  few  reliable  measurements  of  the 
temperature  distribution  in  the  supersonic/hypersonic  flow  regime  or  even  of  the 
turbulence  quantities  themselves.  Existing  turbulence  models  currently  in  use  can 
easily  be  shown  to  be  incompatible  with  experimental  data  as  well  as  inconsistent  in  the 
asymptotic  limit  of  high  Reynolds  number  flows. 

The  objective  of  this  research  was  to  determine  the  asymptotic  structure  of  a 
high-speed  compressible  turbulent  flow  in  the  Umit  of  Itirge  Reynolds  numbers  and 
further  to  develop  simple  turbulence  models  which  can  predict  skin  friction  and  heat 
transfer  coefficients.  The  study  of  He  et  al.  (1992)  reports  the  development  of  a  self- 
consistent  mathematical  asymptotic  structure  for  high-speed  compressible  turbulent 
boundary  layers.  In  this  paper  the  adiabatic  wall  was  considered  and  an  analytical 
expression  was  derived  for  the  recovery  factor.  New  and  simple  turbulence  models  were 
obtained,  and  a  direct  comparison  with  experimental  data  showed  excellent  agreement 
for  both  the  velocity  and  total  enthalpy  data.  Subsequent  work  hcis  incorporated  the 
influence  of  heat  transfer  with  the  new  models  and  also  to  extend  the  theory  into  the 
hypersonic  regimes.  The  theory  described  in  He  et  al.  (1992)  is  valid  for  mainstream 
Mach  numbers  up  to  about  6. 

A. 4  Unsteady  Three-Dimensional  Separation 

Associated  Graduate  Student:  Robert  Puhak 
Status:  Ph.D,  expected  in  December,  1993. 

Recent  adv<mces  in  the  theory  of  unsteady  separation  (Van  Dommelen  & 
Cowley,  1990;  Cowley  et  al.  1991)  have  suggested  that  when  a  three-dimensional 
boundary  layer  undergoes  a  process  of  unsteady  separation  from  a  surface,  a  generic 
structure  develops  in  a  high  Reynolds  number  flow.  It  is  suggested  that  the  eruption 
takes  the  form  of  a  sharply  focussed  tongue  of  fluid  that  moves  rapidly  away  from  the 
surface.  In  three-dimensional  flows,  boundary-layer  sep^u•ation  is  believed  to  initiate 
somewhere  along  the  surface  of  zero  vorticity  (where  both  the  transverse  and 
streamwise  vorticity  components  vanish  simultaneously).  However,  numerical  work  on 
three-dimensional  separated  flows  is  rather  limited.  The  objective  of  this  research  is  to 
calculate  the  evolution  of  some  fundamental  problems  involving  three-dimensional 
unsteady  separation  of  a  surface  boundary-layer.  The  goal  is  to  track  the  evolution  of 
the  different  types  of  recirculations  that  can  develop  and  from  this  determine  how  the 
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layer  evolves  toward  an  eruption. 

A  model  problem  under  consideration  consists  of  a  circular  cylinder  mounted  on 
a  flat  surface  and  a  high  speed  flow  moving  steadily  from  left  to  right.  This  geometry 
represents  a  generic  type  of  configuration  occurring  in  a  large  number  of  situations 
(wing/air  frame  junctures,  turbine  blade  base  flows,  for  example),  <ind  it  is  important  to 
understand  the  dynamics  of  such  regions.  At  present  numerical  solutions  for  the 
developing  flow  on  the  upstream  and  downstream  symmetry  plane  have  been  obtained. 
The  upstream  symmetry  plane  is  of  particular  interest  since  here  the  adverse  pressure 
gradient  due  to  the  cylinder  gives  rise  to  unsteady  separation  in  the  end-wall  boundary- 
layer.  The  present  computations  for  the  developing  boundary  layer  on  the  end  wall  are 
among  the  first  to  be  carried  out  in  three  dimensions  in  Lagrangian  coordinates. 
Calculated  results  show  the  evolution  of  a  complex  separation  singularity  (Degani, 
Puhak  &  Wcilker,  1993). 
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B.  Professor  Varlev 

B.l  Applications  of  Some  Exact  Solutions  to  the  Navier-Stokes  Equations 

Two  families  of  exact  solutions  to  the  Navier-Stokes  equations  (NSE)  have  been 
used  to  analyze  some  technicaJly  important  flows.  The  first  family,  designated  FI 
flows,  is  characterized  by  the  fact  that  the  component  of  fluid  velocity  in  the  direction 
normal  to  some  plane  is  a  function  only  of  time,  t,  and  the  distance,  y,  from  that  plane. 
The  plate  can  rotate  tind  translate  with  time  varying  velocities.  In  addition  to  being 
exact  solutions  to  the  NSE,  the  FI  solutions  are  exact  solutions  to  the  boundciry-layer 
equations  (BLE):  the  terms  that  are  usually  omitted  to  obtain  the  BLE  from  the  NSE 
are  identically  zero.  Well-known  examples  of  parallel  Fl  flows  are  Couette  (1890), 
Rayleigh  (1911)  and  Ekman  (1905)  flows.  Examples  of  non-parallel  flo,vs  are  the 
various  von  Karman  (1921)  swirling  flows,  and  the  Blasius-Hiemenz  (1908,  1911)  flow  at 
a  forward  stagnation  point.  In  these  examples,  either  the  fluid  is  contained  between 
two  parallel  plane  surfaces,  or  occupies  a  semi-infinite  region  bounded  by  a  single  plane 
surface.  For  Couette  and  Rayleigh  flows,  the  bounding  surfaces  translate  in  directions 
parallel  to  their  planes,  but  do  not  rotate;  for  Ekman  flows  the  bounding  surfaces  rotate 
about  the  same  fixed  axis  with  the  same  constant  total  angular  velocity  and  translate  in 
directions  parallel  to  their  planes,  von  Karman  flows  are  generated  when  two  parallel 
plane  surfaces  rotate  with  different  angular  speeds  about  a  common,  fixed  axis  that  is 
normal  to  their  planes.  The  Blasius-Hiemenz  flow  approximates  that  produced  near  a 
line  of  symmetry  when  a  fluid  is  squeezed  between  two  parallel  plates  moving  with 
different  speeds  normal  to  their  planes. 

All  the  viscometric-type  flows  listed  above  axe  produced  when  the  motions  of  the 
bounding  plane  surfaces  have  at  most  two  degrees  of  freedom.  The  Fl  solutions  to  the 
NSE  obtained  in  the  course  of  this  research  may  be  used  to  describe  flows  that  are 
produced  when  the  bounding  surfaces  move  in  any  manner  that  is  compatible  with  the 
remaining  plane,  rigid  and  parallel:  they  c«m  translate  in  directions  pcirallel  and  normal 
to  their  planes  and  rotate  about  different  axes.  The  mathematical  problem  reduces  to 
that  of  solving  a  set  of  PDE  in  the  two  independent  variables  (t,y).  This  can  be  split 
into  two  sets;  the  first  set,  which  is  usually  nonlinear,  describes  the  primary  component 
of  the  flow  and  is  identic«il  to  that  governing  one  of  the  flows  listed  above:  it  uncouples 
from  the  second  set  which,  essentially,  is  linear  with  coefficients  that  tire  determined  by 
the  primtiry  component  of  the  flow.  Two  sub-classes  of  Fl  flows  are  of  special  interest 
because  they  involve  relatively  simple  mathematics.  The  simpler  one  is  characterized 
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by  the  fact  that  the  prims'  ’  component  is  such  that  the  coefficients  in  the  second  set  of 
linear  PDE  are  functions  only  of  t.  These  equations  can  be  transformed  into  two, 
uncoupled,  constant  coefficient  diffusion  equations.  For  the  second  sub-class  of  Fl 
flows,  the  primary  component  of  the  flow  is  steady  and  is  governed  by  a  set  of  nonlinear 
ODE;  the  second  set  of  linear  equations  have  coefficients  that  are  functions  only  of  y. 

The  second  family,  the  F2  flows,  are  generalized  plane  flows  in  the  sense  that,  if 
(xj,  X2,  y)  denote  the  Cartesian  co-ordinates  of  a  point  in  the  fluid  jt  which  the  velocity 
components  are  (uj,  U2,  v),  the  u,  are  independent  of  y.  One  example  is  Poiseuille 
(1984)  pipe  flow,  for  which  the  u^  =  0  tind  v  =  v(X|,  X2);  another  is  the  Burer  (1941) 
vortex  flow,  for  which  v  =  Vjy  and  the  flow  is  radially  symmetric  about  the  y-axis.  As 
another  example,  we  have  analyzed  unsteady,  radially  symmetric,  F2  flows  for  which 
V  =  V|(t,r)y -t- Vglt.r).  These  are  produced  netir  the  outer  surface  of  a  circular  cylinder 
whose  velocity  varies  in  time:  the  surface  may  rotate  with  a  time  varying  angular 
velocity,  translate  with  a  time  varying  velocity  in  the  direction  parallel  to  the  axis  of 
the  cylinder,  and  expand  or  comract  with  a  time  varying  velocity  in  the  radial 
direction.  Away  from  the  cylinder  the  flow  asymptotes  to  that  in  an  unsteady  Burger 
vortex. 

The  main  result  of  this  research  has  been  to  show  how  the  Fl  solutions  of  the 
NSE  may  be  used  to  analyze  some  features  of  the  strong  local  interaction  between  a 
thin  shear  layer,  in  which  the  flow  is  initially  parallel,  and  a  disturbance  in  the  outer 
mainstream.  An  excellent  example  of  such  a  problem  occurs  in  the  study  of  the  strong 
local  atmospheric  disturbance  that  is  produced  in  a  layer  of  air,  where  the  flow  is 
strongly  sheared,  by  the  passage  of  a  tornado.  Here,  the  lower  boundary  of  the  layer, 
representing  the  earth,  is  fixed  and  the  upper  boundary,  representing  the  cloud  base,  is 
rotating  and  translating.  Mainly,  we  have  concentrated  on  the  flow  in  regions  away 
from  material  boundaries  where  it  is  the  fact  that  the  fluid  is  strongly  sheared  that 
controls  the  interaction  with  the  mainstream  flow  and  not  the  presence  of  rigid 
boundaries.  In  this  sense  the  layer  is  a  free  shear  layer,  and  the  governing  equations  can 
be  transformed  into  constant  coefficient  diffusion  equations. 

The  solutions  to  the  NSE  equations  describing  Fl  flows  may  be  used  only  to 
describe  the  local  behavior  of  global  flow  solutions  that  must  be  determined  either 
numerically  or  by  using  approximate  equations.  These  approximate  equations  are 
obtained  by  noting  that  these  global  flows  are  often  well-described  by  solutions  of  the 
BLE  in  which  the  pressure  p(t,  Xj,  X2)  can  be  estimated  a  priori.  Moreover,  away  from 
rigid  boundaries,  even  though  the  flow  is  strongly  sheared,  the  viscous  terms  in  the  BLE 
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are  unimportant  and  can  be  omitted.  We  have  described  a  procedure  that  can  be  used 

to  integrate  the  resulting  hydraulic  flow  equations  (HFE).  This  involves  reformulating 

the  equations  in  terms  of  Lagrangian  variables,  thus  reducing  the  problem  of  solving  a  * 

system  v>t  uoniineM  PDE  to  that  of  solving  a  system  of  nonlinear  ODE.  In  the  special 

case  of  FI  flows,  the  pressure  gradient  terms,  3p/3x,-,  sire  linear  forms  in  the  x-  with 

coefficients  that  depend  on  t.  The  resulting  ODE  are  linear,  and  a  general 

representation  for  the  flow  variables,  which  are  also  exact  solutions  for  the  Euler  * 

equations  (EE),  can  be  found.  These  solutions  are  used  to  show  that  unless  p(t,  Xj,  X2) 

is  estimated  with  care,  the  Fl  solutions  to  both  the  EE  and  the  NSE  develop 

singularities  in  a  finite  time. 

As  an  example  of  our  general  analysis  we  have  described  the  flow  that  is 
produced  at  the  edge  of  a  shear  layer  during  the  passage  of  a  Rankine  vortex  in  the 
mainstream.  Far  from  the  layer  the  axis  of  the  vortex  is  perpendicular  to  the  layer  and 
moves  with  a  time  varying  velocity;  the  strength  of  the  vortex  also  varies  with  time.  ^ 

This  flow  models  that  occurring  in  *^he  center  of  a  tornado  which,  far  above  the  earth,  is 
convected  with  the  ambient  wind-speed.  The  resulting  flow  is  unsteady  and  fully  three 
dimensional. 

The  Fl  and  F2  solutions  of  the  NSE  ase  exact  solutions  to  flow  problems  that  are  $ 

highly  idealized  in  that  they  involve  boundaries  of  infinite  extent.  However,  like  the 
classical  Rayleigh  and  von  Karman  solutions,  sometimes  they  can  be  used  to  provide 
useful  information  about  local  condition  in  more  complicated  and  realistic  global  flows. 

For  example,  Fl  flows  can  be  used  to  approximate  a  real  flow  in  the  vicinity  of  an  axis  » 

of  rotation  or  a  line  of  symmetry,  while  F2  flows  can  be  used  in  the  vicinity  of  some 

plane.  The  Fl  and  F2  solutions  are  similarity  solutions  of  the  NSE  and,  like  other 

similarity  solutions,  must  be  interpreted  with  care:  they  can  be  u.sed  only  when  it  can 

be  argued  that  the  local  flow  is  well  approximated  by  a  viscometric  flow.  I 
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C.  Wave  Induced  Ignition 
C.l  Introduction 

Exothermic  reactions,  with  Arrhenius  kinetics,  depend  exponentially  on  the 
cictivation  energy.  In  many  cases  this  chciracteristic  energy  is  relatively  large  and  rapid 
thermal  heating  occurs.  Over  the  induction  period  prior  to  ignition  the  evolution  of 
small  gas  dynamic  disturbances  is  governed  by  the  Clarke  equation  (Clarke  1981,  see 
§C.2).  Analytical  progress  with  this  equation  is  difficult.  Numerical  solutions  for  semi- 
inhnite  domains  in  which  the  reaction  is  switched  on  by  the  passage  of  a  shock  of  finite 
strength  have  been  obtained  by  Clarke  &  Cant  (1985)  and  by  Jackson  &  Kapila  (1985). 
More  recently,  Blythe  i:  Crighton  (1989)'  developed  an  asymptotic  approach  based  on 
the  Newtonian  limit  »1  where  is  the  frozen  specific  heat  ratio.  Excellent  agreement 
with  the  numerical  calculations  is  found. 

More  complex  problems  arise  when  the  motion  is  generated  in  a  finite  domain  by 
an  initial  (weak)  non-uniform  temp<‘rature  distribution.  The  ignition  stage  of  these 
disturbances  is  a  precursor  to  the  development  of  a  local  hot  spot  into  a  fully  fledged 
thermal  exposition,  and  is  again  governed  by  Clarke’s  equation.  Jackson,  Kapila  &: 
Stewart  (1989)  have  investigated  this  cla.ss  of  problems  by  using  co-ordinate  expansions 
combined  with  computational  techniques.  As  in  the  shock  generated  problem, 
analytical  progress  is  difficult. 

The  research  carried  out  under  this  part  of  the  grant  was  concerned  with 
extending  the  Newtonian  approach  to  examine: 

(1)  non-linear  hot-spot  structure  and  the  reaction  wave  spe«!d  (§C.3); 

(2)  influence  of  upstream  disturbances  on  shock  induced  ignition  (§C.4); 

(3)  spatio-temporal  evolution  of  thermal  irregularities  in  bounded  domains  (§C.5). 

These  problems  have  direct  bearing  on  fundamental  problems  related  to  the 
birth  of  detonation  waves  and  the  onset  of  thermal  explosions  in  confined  vessels. 
Combustion-driven  disturbances  that  propagate  into  non-uniform  regions,  possibly  far 
from  thermodynamic  equilibrium,  also  tirise  in  flow  fields  that  are  characteristic  of 
proposed  hypersonic  transports. 


'Hereafter  referred  to  as  BC. 
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C.2  Iks  Clark?  Euufttiga 

Small  disturbance  theory  at  high  activation  energies  f  “  *  is  governe<l  by  a 
distinguished  limit  in  which  the  disturbance  amplitude  is  0(<).  Here  f  ~  *  is  a 
dimensionless  measure  relative  to  a  background  state  characterized  by  a  (dimensicaless) 
temperature  T  =  1.  Over  the  induction  time  scale  local  solutions  can  be  expwded  in 
the  form 


T  =  To  +  fT‘T,(x,t)  + 


(C.2.1) 


where  Tq  is  either  unity  or  the  initial  temperature  behind  any  reaction- inducing  shock. 
After  suitably  scaling  the  spatial  co-ordinate  x  and  the  time  t,  the  temperature 
perturbation  T |  can  be  shown  to  satisfy  Clarke’s  equation 


V3t^  d^j  \5t^  dx^j 


(C.2.2) 


In  (C.2.2)  the  high  frequency  operator  is  associated  with  the  frozen  wave  speed  7*^^  and 
the  low  frequency  op<‘rator  with  the  isothermal  (Newtonian)  spcH'd  (  =  1  in  current 
variables). 

Although  general  analytical  solutions  of  (C.2.2)  do  not  appear  to  be  possible, 
progress  has  been  made  in  the  Newtonian  limit  7— +1.  To  leading  order  the  wave 
equation  (C.2.2)  degenerates  into 


5x2A5t  I 


(C.2.3) 


which  can  be  integrated.  However,  as  shown  in  BC  for  shock-induced  ignition,  this 
result  is  not  uniformly  valid  near  the  ignition  point  (hot  spot)  where  local  adjustments 
take  place  in  a  region  that  is  exponentially  thin  with  respect  to  7-  1.  Discussion  of  the 
structure  near  the  hot  spot  is  reported  in  §C.3.  Note  that  this  analysis  follows  BC  and 
uses  the  ordered  limit 


1  >7-1 


(C.2.4) 


whereas  the  derivation  of  (C.2.2)  is  dependent  only  on  the  a,ssumption  that  e  <  1. 

Modification  of  the  Newtonian  analysis  given  in  BC  is  cdso  required  when  (C.2.4) 
is  replaced  by 
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■y-l  =0(f)  =  o(l)- 


(C.2.5) 


For  this  limit,  the  Clarke  equation  is  still  valid,  but  the  ignition  time  upstream  of  the 
shock  is  now  comparable  with  that  downstream  of  the  shock.  Consequently,  thermal 
excitation  ahead  of  the  shock  must  be  included  in  the  analysis.  Basic  results,  including 
corrections  to  the  ignition  time,  are  outlined  in  §C.4. 

As  noted  in  the  Introduction  (§C.l),  equation  (C.2.2)  also  governs  the  evolution 
of  initial  spatial  thermal  irregularities  in  bounded  domains.  The  amplification  of  local 
hot  spots  in  such  flows  has  bc-en  examined  by  Jackson,  Kapila  &  Stewart  (1989)  when 
7-l  =  0(l).  An  analytical  approach,  based  on  the  Newtonian  limit  (C.2.4),  is 
described  in  §C.5. 

C.3  Hot  Snot  Structure  and  Wave  Speed 

For  shock  induced  ignition,  the  disturbance  is  initiated  by  a  piston  moving  at 
finite  speed.  Using  expansions  of  the  form  (C.2.1),  the  small  disturbance  equations  can 
be  written  in  characteristic  form  as 

ai(p,±vi)=e^»,  (C.3.1) 

with 

Here  pj  is  the  pressure  perturbation,  Vj  is  the  velocity  perturbation  (relative  to  the 
piston  spt'ed),  ^  is  the  distance  from  the  piston  face,  and  the  time  to  ignition 


7-  =  t,-  -  t. 


(C.3.3) 


where  t,  is  the  ignition  time  on  the  piston  face  and  must  be  determined  as  part  of  the 
solution.  The  exponential  term  on  the  right-hand  side  of  (C.3.1)  represents  the  strong 
influence  of  exothermic  heating.  Similarly,  along  the  particle  paths 


aT,  ap, 


=  il-S) 


with 


(C.3.4) 


6  =  (7-l)/7- 


(C.3.5) 
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These  results  are  based  ou  a  one-step  first  order  reaction  with  Arrhenius  kinetics, 
and  the  corresponding  reduced  rate  law  for  the  perturbation  mass  fraction  yj  is 

=  (C.3.6) 

where  /3  is  a  dimensionless  formation  energy. 

Newtoniai.  expansions  of  the  small  disturbance  equations  correspond  to  the 
subsequent  limit  >0  and,  prior  to  ignition,  have  the  form 

TiU,  r;^)  =  Tio(6r)-h/iT„U,r)-^...  {C.3.7) 

etc.  Substitution  into  (C.3.1)  -  (C.3.6)  shows  that  the  solution  for  the  leading  order 
terms  can  be  written 


Tio=  -/^yio=  -«n(r-(-v-  'O 

Pio  =  ^-^  +  +  +  0  +  +  0  I 

vio  =  «n(T'  + V"  + 


where 


(7-1)M^42 

_27M2-(7-1) 


(C.3.8) 


(C.3.9) 


•% 


» 


»  « 


and  Mq  is  the  initiating  shock  Mach  number.  Determination  of  the  functions  f(7-)  and 
g(r)  requires  implementating  of  suitable  shock  relations  along  the  shock  path 

^  =  V(t.-r).  (C.3.10) 

Solutions  of  the  resulting  functional  equations  can  be  expres.sed  in  the  form 

n  _  oo 

f(r)=  L  (-N)"M1-A"(1  -r)]  (r>0) 

(C.3.11) 

oo 

-^)1  (’•>-V) 

n  =  1 
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with 


A-  1-V 

^"TTV’ 


NT  _  1  —  R 


R  =  V. 


{C.3.12) 


Note  that  1>A>0,  1>N>0. 

Further  details  can  be  found  in  BC,  together  with  determination  of  higher  order 
corrections.  In  particulau,  it  is  shown  that  the  ignition  time 


t^  =  1  +  + ..., 


(C.3.13) 


where 


a  = 


(1+A) 

2  A 


(C.3.14) 


From  (C.3.8)  it  can  be  seen  that  the  solution  has  a  logarithmic  singularity  at 
r  =  0  on  the  piston  face  ^  =  0.  In  fact  (C.3.8)  suggests  that  the  singularity  propagates 
along 


but  this  concerns  events  for  t<0  (t  >t^)  where  the  above  small  disturbance  theory  is 
not  valid. 

The  basic  Newtonian  expansion  (C.3.7)  fails  near  the  hot  spot  ^  =  r  =  0,  and  it  is 
convenient  to  introduce  a  centered  co-ordinate 


» 


X  =  ih- 


(C.3.15) 


From  (C.3.8),  the  limiting  behavior  of  the  outer  temperature  becomes 

Tj~  -toT-(!n(l-)-V-‘x)  +  0(«) 

Newtonian  small  disturbances  theory  also  predicts  that  the  remaining  thermodynamic 
variables,  including  the  density,  have  logarithmic  singularities  on  the  piston  face  at 
T-0.  Near  ^  =  r  =  0  it  is  necessary  to  introduce  (r>0)  a  nonlinear  stretched  time 
through 


(C.3.16) 
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a  =  -  in  c  (t  =  exp{-<T/6))  (C.3.17) 

and  to  seek  local  expansions 

Ti=f^-f(x,  6), 

Pi=^+P(x,<^;^).  }  (C-3-18) 

Vi  =  U(x,  cr;  6), 


etc. 

It  is  not  difficult  to  show  that  the  leading  approximations  to  ^  is  given  by 


-«n[l  +  C(<T)x].  (C.3.19) 

Here  C(<t)  is  defined  implicitly  by 

Ce^  =  V ~  ^  exp  [V ~  ^ -a],  (C.3.20) 


and  the  derivation  of  (C.3.20)  requires  extending  the  analysis  to  second  order  terms  (see 
BC).  Related  results  for  the  remaining  dependent  variables  can  also  be  obtained.  In 
particular,  the  leading  order  velocity 


c 

l-C^ 


in 


l  +  X 

.1+CxJ’ 


(C.3.21) 


which,  as  required,  vanishes  on  X  =  0-  Further,  the  logarithmic  singularity  in  the  outer 
density  is  replaced  by  a  bounded  limit  as  <r— ►oo(t— >0),  which  is  consistent  with  the 
numericed  cadculations  of  Jackson  &  Kapila  and  Cl£irke  &  Caint. 

The  representation  (C.3.18)  is  not  valid  across  r  =  0,  i.e.  in  r  <  0  (see  (C.3.17)). 
To  extend  the  solution  into  r  <  0  it  is  necess£u^y  to  introduce 


t  =  tI(  = 


i 

X' 


(C.3.22) 


so  that  in  T  >  0  (r  >  0) 


» 


» 
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<r  =  s  -  ^  to  r. 


(C.3.23) 


Either  directly  from  the  governing  equation,  or  by  re-expressing  (C.3.19)  -  (C.3.21)  etc. 
in  terms  of  s  and  r,  it  can  be  shown  that 


T,  =  I  -  to  l,  .  C(s)l  +  ^  I  {r  -  to  (1  +  r)}  +  r,(,) 


+  .  . 


(C.3.24) 


with  corresponding  results  for  the  other  dependent  variables.  The  expression  (C.3.24)  is 
valid  in  r  <  0(t  <  0).  In  general,  this  form  of  the  solution  is  not  uniformly  valid  as  r— »oo 
(i.e.  on  the  piston  face  with  r<0),  unlike  the  formulation  (C.3.18),  which  is  valid  on 
^  =  0  for  T  <  0.  The  representations  in  terms  of  (s,  r)  and  (ct,  x)  do,  however,  overlap 
for  r  =  0(1). 

From  (C.3.24),  it  appears  that  the  ignition  path  is  defined  by  the  logarithmic 
singularity  on 


r  =  r,-=  -C(s), 


(C.3.25) 


cJthough  the  term  O(^)  has  a  stronger  algebraic  singularity  along  this  path.  Results 
obtained  by  Friedman  &  Herrero  (1989)  imply  that  the  Clarke  equation  can  at  most 
have  a  logarithmic  singularity  at  ignition.  This  suggests  using  a  PLK  approach  to  define 
ignition  through 


r.=  -C(s) 


(C.3.26) 


with 


s  =  s  -I-  ^Sj(s)  +  .  .  . 

S|(s)  is  now  chosen  to  eliminate  the  algebraic  singularity  which  requires 

,,=  .C-to(l-C).ii±^. 


(C.3.27) 


(C.3.28) 


Using  this  result,  and  (C.3.26),  it  can  be  established  that  the  reaction  wave  speed 


li-'C 


rr^Mi^l 
^  \  C  l  +  C  J  + 


(C.3.29) 
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where  now  C  =  C(s).  Note  that  the  wave  speed  becomes  unbounded  on  the  piston  face 
(^  =  0  or  s  =  oo,  implying  C  =  0);  this  is  consistent  with  numerical  prediction  for  the 
Clarke  equation. 

Internal  to  the  present  region,  which  is  exponentially  thin  with  respect  to  6,  the 
complete  structure  for  the  full  equations  includes  a  fuel  consumption  zone  that  is 
exponentially  thin  in  e.  Corresponding  extensions  for  the  full  wave  structure  tire 
planned. 

The  result  (C.3.28)  also  fails  near  C  =  1  where  the  wave  speed  is  equal  to  the 
speed  of  the  outgoing  char<icteristic  leaving  the  hot  spot.  Note  also  that  the 
temperature  solution  (C.3.24)  has  a  logarithmic  singularity  on  this  characteristic 
(r  =  - 1),  but  the  present  inner  solution  is  valid  only  in  r  >  -  C(s)  with  C  <  1. 

Resolution  of  these  difficulties  requires  an  analysis  of  the  full  equations,  including 
fuel  consumption,  2uid  is  an  essential  step  in  providing  a  complete  description  of  the 
ultimate  evolution  of  any  detonation  wave.  Future  research  will  examine  this  structure. 

C.4  Upstream  Heat  Release 

If  the  ordered  limit  (C.2.4)  is  replaced  by 

1>7-1=0(€),  (C.4.1) 

then  modifications  must  be  made  to  the  analysis  described  in  §C.3.  Using  the  expansion 
(C.2.1),  with 

y  =  l+cT2yj(^,  r)  +  .  .  .  ,  (C.4.2) 


the  full  rate  law 


» 


» 


» 


dry  =  ^  ^ 


y  exp 


(C.4.3) 

» 


reduces  to  (C.3.6)  downstream  of  the  shock.  In  generail,  exothermic  heating  will  also 
occur  upstream  of  the  shock.  Suitable  local  expansions  are 


» 


•  •  •  • 


•  •  • 
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Tj  -  1  +<Tjg^(t)  + .  .  . 


yi  =  1  +fyioo(t)  +  -  •  ■ 


(C.4.4) 


with  pjgp  =  Tjoj,  and  =  0.  It  is  readily  established  that 


i2e  _  1  p^ioo 


3t  tgQ 


(C.4.5) 


where 


‘oo  =  exp  ? 


'(¥)1 


(C.4.6) 


When  7  - 1  =  0(f),  the  standard  shock  relations  imply  that  T(,-l=0(f),  so  that 
t^=:0(l),  i.e.  ignition  times  upstre<im  of  the  shock  are  now  comparable  with  those 
downstream  of  the  shock. 

For  convenience,  it  is  supposed  that  the  shock  is  initiated  at  t  =  t,,  and 
propagates  into  the  region  x  >  0.  The  upstream  reaction  is  assumed  to  start  at  t  =  0 
with  an  initial  state 

T  =  p  =  p  =  y  =  l,  u  =  0. 


Consequently,  from(C.4.5) 


(C.4.7) 


Obviously,  for  the  current  analysis  to  be  appropriate,  to^  >  to- 


The  solution  downstream  of  the  shock  is  still  governed  by  the  Clarke  equation, 
but  conditions  on  the  shock  path 


^  =  V(t--r),  (t,-  =  t,--to) 


(C.4.8) 


must  be  modified.  In  particular,  the  relationship  between  the  temperature  and  velocity 
perturbations  is 


•  • 


(C.4.9) 


T^i  -<5Qv,i  +  QooTi^, 


where 


2(7M^-H)V  _ 

{M2  +  1)((7-1)M2  +  2))’ 


(C.4.10) 


Qoo  ~  1  ^Qooi 


(C.4.11) 


and 


27  (Mg  -  1  )(7(7  + 1  7(7  -  1)M^  -h  7(7  +  ml  +  2) 

(7  +  1)2  m2[(7-1)m2  +  2] 


(C.4.12) 


Similarly, 


Pfil  =  +l^ooTioo' 


(C.4.13) 


with 


» 


(37  +  1)M2-(7-1) 
*  (7  +  1)^oT2(M2  +  1))' 


(C.4.14) 


Here  the  coefficients  denoted  by  a  subscript  c?o  are  associated  with  the  upstream 
behavior. 

From  these  results,  the  leading  order  solution  for  the  temperature  perturbation  is 
found  to  be 

Ti= -to  [r  +  (l-t-l)V-l^].  (C.4.15) 

In  addition,  relative  to  shock  initiation,  the  ignition  time  is 

t,-to  =  l-^.  (C.4.16) 

*■00 

Corresponding  results  for  the  velocity  and  pressure  perturbation  have  been  determined 


» 


» 


» 


•  •  •  • 
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and  can  be  written  down  in  a  form  similar  to  (C.3.8)  and  (C.3.11).  Similarly, 
extensions  of  the  analysis  to  describe  the  hot  spot  structure  have  been  carried  out  using 
the  approach  described  in  §C.3.  Calculations  of  higher  order  corrections  to  the 
temperature  distribution  and  the  ignition  time  are  partially  complete. 


C.5  Ignition  frcnn  a  Spatially  Distributed  Source 

(a)  Basic  Theory 

For  a  spatially  uniform  initial  state,  ignition  should  ideally  occur  in  a  spatitdly 
homogeneous  manner.  Weak  irregularities  in  the  background  thermal  field,  however, 
can  lead  to  local  hot  spots  that  become  energetic  centers  associated  with  the 
development  of  a  thermal  explosion.  Problems  of  this  type  have  been  considered  for 
bounded  domains  (0  <  x  <  ^)  by  Jackson,  Kapila  &  Stewart  (1989)^.  Over  the  induction 
phase,  prior  to  ignition,  it  can  be  shown  that  the  disturbance  is  still  governed  by 
Clarke’s  equation. 

In  this  section  the  variables  adopted  will  follow  those  introduced  by  JKS. 
Relative  to  (C.2.2),  the  coordinates  (x,t)  are  now  replaced  by  (7^^^x,  ^t)  and  the  Clarke 
equation  becomes 


Similarly,  the  results  (C.3.1)  and  (C.3.4)  can  be  re-written  ais 


(C.5.1) 


(C.5.2) 


5Ti  ^ 
dt  dt 


(C.5.3) 


with,  as  before,  6  =  It  is  assumed  that  there  is  an  initially  thermal  non¬ 

uniformity  of  size  e  so  that 


Ti(x,O)  =  0(x/£) 

with  i  (C.5.4) 

Pl(x,0)  =  Ui(x,0)  =  0. 


^Hereafter  referred  to  as  JKS. 


» 


» 


» 
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In  addition,  the  velocity  must  satisfy 


uj(0,t)  =U2(£,t)  =  0. 


(C.5.5) 


For  convenience,  the  domain  0  <  x  <i  is  mapped  onto  the  unit  interval  by 
introducing 


(x,t)=>£(x,t), 


(C.5.6) 


together  with 


Tj=>T2-£n^,  pj=>Pi+£n£ 


Pl=>Pl,  Ui=>U2. 


(C.5.7) 


Under  this  mapping,  the  basic  equations  are  invariant.  Boundary  and  initial  data  are 
also  unchanged  save  that 


Tj(x,0)  =  <),(x)  +  «n  £  =  ^(x;e). 


(C.5.8) 


Newtonian  expansions  of  the  form 


Tj  =  TjQ(x,t) +  <5Tj2(x,t)  +  ... 


(C.5.9) 


etc.,  are  used  to  generate  the  small  amplitude  pre-ignition  solution.  As  in  §3,  it  is 
iissumed  that  1  >  ^  >  e.  The  basic  solutions  for  the  temperature  is  easily  shown  to  be 


Ti^=  -£n[te(x)-tj. 


(C.5.10) 


where  the  appcirent  ignition  locus 


te(x)  =exp{  -«(x;«)}. 


(C.5.11) 


Further  analysis  then  gives 
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4Pio  =  +  I  +  (x,«)  -  I  -  (X,^i) 

4Uj^  =  f(Q)-f(;3)  +  I  +  (x,o)  +  I_(x,  li) 
where 


ds 

t^(s)  -  (r±s) 


cuid 

a  =  t-x,  (3  =  t  +  x. 


}■ 


(C.5.12) 


(C.5.13) 


(C.5.14) 


The  function  f(t)  can  lx‘  fouiul  hy  application  of  the  initial  and  iKnindary 
conditions.  It  follows  that  for  2n  >  t  >  2n  -  I, 


f(t)  =  .I^t  -  (2ni  -  1 ))  -  I  +  (2n  -  t,  t  -  2n), 

m  =  1 

and  for  2n  +  1  >  t  >  2n, 

f(t)  =  J  j(t-(2rn-l))  +  I_(t-2n,  t-2n). 


(C.5.15) 


(C.5.’6) 


m  =  1 


where 


J(t)  =4[‘  - (is. 

Jo(Us)-t)2-  l-s)2 


(C.5.17) 


Note  that  for  -  1  <  t  <  0 


f(t)=  -I  +  (-t,t).  (C.5.18) 

This  definition  of  f(t)  can  be  n.sed  to  provide  a  complete  representation  for  the  leading 
order  solution. 

As  in  the  shock  driven  ca.se,  it  is  possible  to  obtain  a  higher  order  correction  for 
the  temp)erature  distribution.  This  result  al.so  yields  an  important  correction  to  the 
ignition  time  (see  §C.5b).  From  (C.5.3),  and  the  expansion  (C.5.9), 


■ 


4 


» 


» 


» 


t 


I 
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Using  (C.5.10),  and  Tj(x,0)  =  pjj,(x,0)  =  0,  leads  to 
Tii(x,t)  =  Pi„(x,t)+^^^  Plo(x,s)ds. 


(C.5.19) 


The  solution  defined  by  (C.5.10)  -  (C.5.19)  fails  near  the  hot  spot  on  x  =  0.  Local 
solutions  of  the  type  discussed  in  §3  for  the  shock  generated  case  c^ul  again  be  obtained 
in  zones  that  are  exponentially  close  to  the  hot  spot.  Details  are  given  in  Crighton  & 
Blythe  (1993).  As  in  §3,  the  analysis  described  above  holds  only  for  t  <  t,-,  where  t,-  is 
the  minimum  ignition  time  on  the  left  end  wall.  Evolution  of  the  reaction  wave  away 
from  the  hot  spot  is  not  described  by  the  Clarke  equation,  and  (C.5.11)  does  not  provide 
an  overall  description  of  the  ignition  path.  Results  for  the  minimum  ignition  time  on 
X  =  0  are  presented  b..iow  in  §C.5(b). 

(b)  Ignition  times 

Initial  thermal  disturbances  considered  by  JKS  were  of  two  types: 

Type  B  is  characterized  by  a  finite  heat  flux  temperature  gradient  at  the  left  hand 
boundary,  where  the  temperature  has  its  maximum  value.  In  the  present  rotation 


^(x;^)  =  a(l  - x)  +  ^n  6, 


(C.5.20) 


where  a  is  a  positive  constant.  From  (C.5.11)  it  follows  that 


te  =  iexp(-a(l-x)). 


(C.5.21) 


with  a  minimum  ignition  time 


t  -ie~“ 
''60  ~ 


(C.5.22) 


on  X  =  0. 


Type  I  is  characterized  by  zero  heat  flux  at  the  left  hand  boundau-y  where  again  the 
temperature  has  its  maximum  value.  For  this  case  JKS  take 
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^(x;^)  =  a{l  -  x^)  +  in  t, 


(C.5.23) 


for  which 


te(x)  =  Jexp  (-a(l-x2)] 


(C.5.24) 


and  again  is  the  minimum  ignition  time  on  x  =  0. 

Both  Ccises  correspond  to  initiating  hot  spots  at  x  =  r  =  0.  As  in  §3,  the  ignition 
time  on  X  =  0  can  be  found  from  the  expansion 


+  +  •  •  • 


(C.5.25) 


Determination  of  t^j  requires  an  analysis  of  the  behavior  of  (C.5.19)  near  the  hot  sp)ot. 
Baised  on  the  results  of  Friedman  A:  Herrera,  it  is  required  that  the  local  singular 
behavior  is,  at  most,  logarithmic.  After  some  analysis  (Crighton  &  Blythe,  1993)  it  can 
be  established  that 

tfii  =  -  I  tjo  (0,s)  ds.  (C.5.27) 

If  the  constants  i  and  a  ^lre  such  that  tg^  <  (0  <  t  <  1) 

For  type  B  profiles  JKS  give,  for  a  =  0.5  jmd  i  =  0.8, 
t  ■  =  0.447, 

while  the  present  analyticad  approach  yields 
t-  =  0.461  +0(^2) 

Both  of  these  results  are  based  on  7  =  7/5  so  that  6  =  2/7.  In  view  of  the  size  of  6,  the 
agreement  between  (C.5.29)  and  (C.5.30)  is  remarkably  good.  Comparison  between  the 
predicted  temperature  distribution  and  JKS  shows  even  better  agreement.  Details  are 
provided  in  Crighton  &  Blythe  (1993).  CaJculations  for  type  I  profiles  are  currently  in 
progress. 


(C.5.29) 


(C.5.30) 
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